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Recent studies introduced biased (degree- dependent) edge percolation as a model for failures in real-life 
systems. In this work, such process is applied to networks consisting of two types of nodes with edges 
running only between nodes of unlike type. Such bipartite graphs appear in many social networks, for 
instance in affiliation networks and in sexual contact networks in which both types of nodes show the 
scale-free characteristic for the degree distribution. During the depreciation process, an edge between 
nodes with degrees k and q is retained with probability proportional to (kq)~ a , where a is positive so that 
links between hubs are more prone to failure. The removal process is studied analytically by introducing 
a generating functions theory. We deduce exact self-consistent equations describing the system at a 
macroscopic level and discuss the percolation transition. Critical exponents are obtained by exploiting 
the Fortuin-Kasteleyn construction which provides a link between our model and a limit of the Potts 
model. 



1. INTRODUCTION 

Scale-free self-similar structures occupy a prominent place among Nihat Berker's research tools. He 
pioneered an interesting class of hierarchical sets, commonly called "Berker lattices" [112] on which a 
family of renormalization group transformations is satisfied exactly. Further, the geometrical critical 
phenomenon of percolation [3] and, recently, the cooperative behavior of interacting degrees of freedom 
on random scale-free networks [4] has enjoyed his thorough attention. In this paper, we discuss a subtle 
variant of a percolation process on scale-free graphs of bipartite structure. 

In recent decades, a detailed analysis of many real-life structures revealed the existence of a common 
property. Studies of social, technological and biological networks have indicated that the probability 
distribution P{k) of the degree k of a node, i.e., the number of links attached to the node, follows 
a decreasing power law P(k) oc fc -7 for large values of k. Notable examples of such so-called scale-free 
networks are the network of co-authors of scientific articles, protein interaction networks, power grids and 
transportation networks. An accurate determination of the topological exponent (or "degree exponent") 
7 in real-life networks is far from easy, but most studies indicate a value of 7 in the range 2 < 7 < 3.5 [5]. 

Most properties of scale-free networks are determined by the existence of a small number of nodes 
which have a large number of links, commonly called hubs. An important consequence of the existence 
of these hubs is the resilience of the network against random removal of links or nodes. Such destruction 
processes are studied in the percolation problem, which mostly focusses on the behavior of the largest 
cluster in the network if links or nodes are removed. When all links and nodes are present in a (infinite) 
network, a finite fraction of nodes will belong to the largest or "giant" cluster. However, when more and 
more links are removed, a point will be reached at which the fraction of nodes contained in the giant 
cluster vanishes. This point, which shows all characteristics of a second-order phase transition, is called 
the percolation transition. If the breakdown of the giant cluster occurs when a finite fraction of links is 
still present in the network as a whole, the network is called fragile. However, if after the removal of an 
arbitrary fraction of links the giant cluster is still present, the network is called robust. 
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Studies of percolation on scale- free networks have revealed that the behavior of the giant cluster depends 
on the topological exponent 7 |6I7) . If 7 > 3, the giant cluster breaks down at a finite fraction of links, 
thus the network is fragile. However, when 7 < 3, the network is robust. Most real-life networks are thus 
robust against random removal of links. Among the numerous applications, the percolation problem has 
for instance been studied as a model for failures in the Internet structure and in terrorist networks [6 , 
or as a model for virus spreading on social networks [5]. 

Not all networks can be described using a single type of vertex. For instance, the network of hetero- 
sexual contacts can best be modeled by two types of vertices, men and women, where links always run 
between nodes of unlike type. Such networks are called bipartite networks. They appear in many social 
structures, for instance in affiliation networks, i.e., networks of individuals joined by common membership 
of communities. The two types of vertices then represent the individuals and the groups, while the links 
between them indicate group membership [5]. One could for instance model the research interests of 
physicists: the first type of nodes then represents the researchers, while the second type of nodes consists 
of the research topics (for instance classified according to the PACS). Most studies concerning bipartite 
networks focussed on the architecture and the building process of networks. This is for example the 
case for sexual-contact networks [S], for listening habits and music genres [10) or for general affiliation 
networks [11112) . Among these studies, it was revealed that in many bipartite networks, the degree 
distribution of both types of nodes shows the scale-free characteristic, for instance in the hetero-sexual 
contacts network [13] . Some studies also discuss the behavior of the network under random removal of 
its links or nodes [14] , 

In recent work, random link removal on monopartite graphs was extended to biased percolation, in 
which links are removed according to the degree of their neighboring nodes [15] . More specifically, to 
each link is attached a weight 



where k and q are the neighboring node degrees and a is the bias exponent which we take to be positive 
(and < 1) for our present purposes [31]. During the depreciation process, the probability to retain a link 
is proportional to the weight of the link. Since a > 0, links between hubs are more prone to failure. 
Therefore, the percolation process is called centrally biased. Such a process is inherent to many social 
networks, where friendships between people with many acquaintances are expected to be weaker and 
last less long than friendships involving people with few connections [TB]. In previous articles, Refs. [TS] 
and [T7] , we showed that the critical behavior of such a biased percolation on a scale-free network can be 
mapped onto a random removal process on another scale-free network. More precisely, biased percolation 
with bias exponent a > on a network with topological exponent 7 belongs to the same universality 
class as random percolation (a = 0) on a network with topological exponent 7, given by 



Therefore, the network is, in the macroscopic limit, robust against biased removal of links as long as 
2 < 7 < 3. If 7 > 3, the network is fragile for the central bias process. Note that 7 > 7. 

The paper at hand discusses biased percolation on bipartite scale-free networks. In Sect. 2 we introduce 
bipartite scale-free networks and present a basic description of a link removal process on those networks. 
The third section introduces the theory of generating functions, which we extend so as to describe bipartite 
networks. The percolation threshold and the scaling of the critical point are extracted from the theory. 
In Sect. 4, the equivalence between percolation and a limit of the Potts model is elaborated. Using 
this equivalence, we construct a finite-size scaling theory which enables the calculation of the critical 
exponents for the percolation transition. In Sect. 5 we present our conclusions. 

2. THE MODEL 

We start from a bipartite network with N nodes, divided in Na nodes of type A and Nb nodes of type 
B, each type with its own degree distribution P A (k) and P B {k). Both degree distributions are assumed 
to follow a decreasing power law and thus are scale free. More precisely, 



w kq = (kq) 



—a 



(i) 



P\k) = Cifc-T* 



(3) 
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for values of k between the minimal and maximal degrees, rrii and Ki, respectively, for i = A,B. In 
Eq. @, Ci is a normalization constant and the exponent 7^ is assumed to be larger than two to ensure 
a finite mean degree. Links exist only between nodes of unlike type. Moreover, we assume that no 
degree correlations or mixing patterns between the nodes, as for instance assortative mixing, occur. The 
probability P^{k) that a randomly chosen edge emerging from a node of type A leads to a type-B node 
with degree k thus only depends on the degree distribution P B {k). More precisely, 

P n(k) = ^^, (4) 

where (■),; denotes the average over the nodes of type i, obtained by using the degree distribution P 2 (k). 
Note that the total number of links attached to the nodes of type A equals the total number of links 
attached to type-i? nodes, i.e. (k)ANA = (QbNb- 

In our percolation process, a fraction / of the links is removed in a single sweep; we call this the 
simultaneous process. An edge between nodes with degrees k and q is retained with probability 

f Wkq , . 

and is removed with probability 1 — p^ q . The weight Wk q is defined in Eq. ([1]) and (w) e denotes the 
average weight of an edge, 

{k)A{k)B 

Note that the random link removal process is recovered if a = 0. The positive bias exponent a should 
be smaller than one in order for the depreciated network to be scale free. Moreover, the depreciation 
process is only well-defined if pk q < 1 for all possible values of k and q. Therefore, our percolation can 
only be used correctly for values of / for which |35j 

f < fu = (w) e (m A m B ) a - (7) 

We now introduce the marginal distribution p\(f) as the mean probability that an edge connected to 
a node of type i with degree k is present when a fraction / of links is reincluded in the network. One 
finds 

rfifl = <^-. M 

Note that the marginal distribution p£ does not depend on the distribution of the nodes of type B and 
vice versa. Using the marginal distribution, the degree distribution P 1 (k) and nearest- neighbor degree 
distribution P l n {k) of nodes of type i in the diluted network can be deduced. One arrives at 

P\k) = f>W I W) ¥ (l-Pl) fe ^ (9a) 

k=k 

The importance of these expressions will become clear when introducing the generating functions in the 
next section. Moreover, it can be shown that 

fp kq = PkPg, (10) 
and thus, according to the arguments in Ref. [17] , the diluted network is still uncorrelated. 
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3. GENERATING FUNCTIONS 

Percolation is often studied using the generating functions approach. By this method, self-consistent 
equations for the size of the giant cluster can be obtained easily. Moreover, the method is exact if the 
diluted network is uncorrelated and if loops in the finite clusters can be ignored, which is justified for 
scale- free networks [IB]. We first briefly introduce the general scheme, closely following the approach of 
Newman [19] . Then, the generating functions method is exploited for the case of biased percolation on 
bipartite networks. 

3.1. Self-consistent equations 

Generating functions are used in a plethora of problems concerning series. The generating function of 
a sequence is the power series which has as coefficients the elements of the sequence [3D]. In percolation 
problems, generating functions that generate the probability distributions characterizing the network arc 
widely used |19I21) . The most important functions are those that generate the degree distribution and 
those that generate the nearest-neighbor distribution. To study the percolation problem, both generating 
functions will be defined in the diluted network, for both vertex types. 

The generating function for the degree distribution P l (k) is defined as 

Ki 

FZ(h) = '£p i (k)e- 
I=i 

while the distribution of the residual edges in the diluted network is generated by 

Ki 

Fi(h)=J2P l n(k)e- h(k - 1 \ 

k=l 

where i — A,B. Substituting Eqs. and ((!]), we obtain 

Ki 

F*(h) = P l (k)(l-pl + e- h pl) k , 

k — in i 

Ft(h) = £ - p\ + e- h P %) k -\ 

k—rrii 

with i = A,B. For our interest, the most relevant generating functions for the percolation problem are 
the ones associated with the probability distribution of the size of the finite clusters, since these quantities 
can be related to the size of the giant cluster. Let Hq generate the probability that a randomly chosen 
node of type i belongs to a cluster of a given finite size. Furthermore, we introduce H\ as the generating 
function for the probability that upon following a randomly chosen edge emerging from a node i towards 
the endnode of type j, a cluster of given (finite) size is reached. If the finite clusters can be treated as 
trees and the diluted network is uncorrelated, these generating functions satisfy coupled self-consistency 
equations, analogous to those derived for monopartite graphs in Ref. |19j . For bipartite graphs, we obtain 



(11) 



(12) 



(13a) 
(13b) 



Hf(h) = e- h F*[H*{h)l (14a) 

fff(ft) = e- h Ff[Hf{h)], (14b) 

Hf(h) = e- h F A [Hf(h)}, (14c) 

H*(h) = e- h F Q B [H?(h)}. (14d) 

Here the function F{[H\(h)] denotes the function F{ wherein e~ h is replaced by H\(h) with i = A, B. The 
percolation threshold can now be derived with the aid of these generating functions and self-consistent 
relations. 

3.2. Percolation threshold 

The percolation threshold is most easily studied by introducing the average cluster size in the diluted 
network, i.e., the average fraction of nodes of type % in a cluster, denoted by S l . Using the properties of 
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generating functions, S z can be related to Hq .20 : 

5* = -flS(0), (15) 

where the dot represents differentiation with respect to h. An expression for the average cluster size S, 
i.e., the average fraction of nodes of type A and B in a finite cluster, is then easily found: 

S = -\ E ^o(O). (16) 

i=A,B 

The average cluster size S in the diluted network can be further worked out by differentiating Eqs. (|13al) 
and (|13bj) with respect to h: 

s = 1 f (k) A (l-FF(0)) + (kMl-Ff(0)) 
2 l-Ff(Q)Ff(Q) 

Hence the average cluster size diverges when 

l = i^(0)Ff(0). (18) 

The percolation criterion (|18l) is the extension for bipartite graphs of the Molloy-Reed criterion, which 
provides a condition for the existence of a giant cluster in a network 6 22]. A similar expression for 
the percolation threshold was already found in studies concerning percolation on general multipartite 
networks [13] . 

Using Eqs. ([8]) and (|13bp an explicit criterion for the critical fraction f c can be found, 



fc = >Jtff?, (19) 

where /* denotes the critical fraction of a monopartitc graph consisting only of nodes of type i with 
degree distribution P l (k), 

p — )Z Li (20) 

Jc <A>i«fc 2 - ito ><-<* 1 - 2a }0' 

Note that Eq. ([T9]l reduces to the criterion for monopartite graphs if ja = 7b, as it should be. Moreover, 
the critical fraction of a bipartite graph can easily be found if the critical fraction of the monopartite 
graphs consisting of only nodes of type A and B is known. As a consequence, also the resilience against 
biased failures of the bipartite network in the macroscopic limit is known completely if the behavior of 
the different subgraphs is known. If one of the two subgraphs is robust against biased percolation, also 
the bipartite network will be robust. The robustness criterion is therefore: 

min(7 A ,7 B ) < 3, (21) 

where j i is defined as in Eq. ([2]) . This is illustrated in Fig. [T] where the robust regimes are indicated by 
the blank (undotted) regions. Only if both 7^ and 7 B are larger than 3, the network is fragile, indicated 
by the regimes dotted in red on Fig. 1. 

For the robust regimes, we can quite easily determine the scaling relation of f c as a function of the 
network size, f c oc iV s , by explicitly evaluating the expectation values in Eq. (I2TJ1) . Replacing the sums 
over the degree distibution by integrals, it can be seen that the moments (k)i and (k 1 ~°')i never diverge 
if central bias is applied to a network with 7$ > 2. The moments {k 2 ~ 2a )i and (k 1 ~ 2a ) i , on the other 
hand, may diverge, but the former will always grow faster than the latter. The behavior of the critical 
fraction as a function of the maximal degrees, K{ for i = A,B, thus stems from the first term, which for 
2 < 7i < 3 diverges as ^- 2a -~i\ Since Ki cx TV- 1 ^ 7 ' -1 ^ [6], the scaling of the critical point as a function 
of the network size can be determined. If both 74 and j B are smaller than 3 (regimes IV and IV on 
Fig. [IJ, <^ is given by 
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Note that this reduces to the result for monopartite graphs, determined in Ref. [15], if ja — 1b- If 
2 < ja < 3 < 7s (regime III on Fig. 1), the exponent is the same as the exponent for a robust 
monopartite graph of type A: 

^- 3 ~^ (23) 



Note that the exponent in the regime III' (2 < jb < 3 < 7a) can be found easily by interchanging 7^4 
and 7s in Eq. 

In sum, we have now provided an extension of the Molloy-Reed criterion for the critical threshold of 
biased percolation on bipartite scale-free networks. Moreover, the scaling of the critical fraction as a 
function of the network size was determined for robust graphs. Our results show that the critical fraction 
of a bipartite network is governed by the critical percolation behavior of the scale-free graphs consisting 
of only one type of nodes. 

4. CRITICAL BEHAVIOR 

In the following section, we introduce and calculate the critical exponents of the percolation transition. 
The first part briefly introduces the Fortuin-Kasteleyn construction, which provides a link between per- 
colation and a limit of the Potts model. Exploiting this link with the Potts model, we can define critical 
exponents for the percolation problem. In the second part of the section, a finite-size scaling theory is 
constructed, in order to calculate the critical exponents in the third part. 
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4.1. Fortuin-Kasteleyn construction 

There exists an exact equivalence between random edge percolation and the q — > 1 limit of the g-state 
Potts model, originally worked out by Fortuin and Kasteleyn in Ref. [23]. Moreover, their proof can easily 
be generalized to incorporate edge- dependent coupling constants in the Potts model and edge-dependent 
removal in the percolation model, respectively. The Fortuin-Kasteleyn construction states that the free 
energy of the q — » 1 limit of the g-state Potts model is the same as the "free energy" of the percolation 
problem [24] . The latter is defined as the generating function of the cluster-size distribution n s , 

nf,h) = {T i n a e- h '\. (24) 

Here the average is performed over all networks in which the probability to retain the edge between nodes 
with degrees k and q is pt q - We can immediately identify the probability Voo for a node of any type to 
be in the infinite cluster and the average cluster size 5, as a function of the fraction of removed links: 



d 2 F 



(25a) 

h=0 



S(J) = 



(25b) 

h=0 



dh 2 

Since we are interested in the behavior close to criticality, we introduce the parameter 

e = f - U (26) 

The usual critical exponents a, j p , (3 and 6 can now be defined for the percolation problem using the 
analogy with the Potts model: 

T(e,0) ~ e 2 - a , (27a) 

Poo(e) ~ e~ 7p , (27b) 

5(e) ~ e?, (27c) 

h 1/s - 1. (27d) 



~dh 



e=0 

Relations among these exponents can be found using a scaling theory. 
4.2. Scaling Theory 

In the following section, we introduce finite-size scaling in order to find critical exponents near the 
percolation transition. In order to solve the scaling relation, we use a Landau-like theory which we derive 
from the exact self-consistent relations, Eqs. (fl"4]) . We closely follow the approach presented in Ref. [T7] . 
which is based on Refs. [25|26l27j . 

Our scaling theory consists of two basic scaling relations. The free energy J 7 of a large but finite 
network with ./V nodes close to criticality can be written in the general form 28 : 

J-(e, h) = N^F (tN 1 /^ , hN 1 '^ , (28) 

where F is a well-behaved function. Close to the critical point, the free energy then scales as 

J"(e,0) oc e v \ (29a) 
F(0,h) oc h Vh . (29b) 

As a second ansatz, the scaling of the cluster size distribution n s can in the thermodynamic limit be 
written as 



n.(e) = s- T G(es' T ). 



(30) 
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Using the scaling forms of Eqs. (128)) and (1301) . standard techniques provide us with exponent relations by 

which all critical exponents can be related to and v t . One arrives at |28j : 

(3 = v e (l-v^), (31a) 

7p = ^(2^-1), (31b) 

a = 2-2/3- 7p , (31c) 

a = (/3 + 7p r\ (31d) 

r = 2 + /3(/3 + 7p )- 1 , (31e) 

5 = (P + 7p)/P- (31f) 

The problem we are left with is to find the exponents Vh and v e for percolation on bipartite scale-free 
networks. 

4.3. Critical exponents 

Using the self-consistent equations, the exponents Vh and v e can be determined in an exact way. We 
introduce the parameters 

4i(e,h) = l-Hi(e,h), (32) 

where i — A, B. As we are merely interested in the behavior near the transition where £< 1, /i< 1 and 
V^i <C 1, we can expand Eqs. (|14al) and (Il4b|) . 

4.3.1. Fragile networks 

We first discuss all regimes in which both j A and j B are larger than three, i.e., the regimes in which 
the network is fragile and has a finite f c . An expansion of the self-consistent equations, Eqs. (1141) . yields 

h-^B = -c?(f c + e)yj A + c4yj 2 A + ... + cf(<p A y-- 2 + ..., (33a) 
h-ipA = -cf(/ c + e)^ B + cfv4 + --- + cf(V'Br B - 2 + ..., (33b) 

wherein all c-constants are positive and according to the percolation criterion f^c^cf = 1. The existence 
of the correction terms, c l s (ip A )' 1i ~' 2 , can be checked numerically. Eqs. (|33|) are derivable by minimization 
of the free energy 

f - *,* B + e (--*. - + M + 4(A + , . . , (34 ) 



with respect to the parameters ipi. Using the equations of state, Eqs. (|33|) . we may write at the saddle 
point 

i— a,b 



z = ?(-?-T + ^f^ + ^'" 1 )- <»> 



If 7^4 and "f B are both larger than 4, i.e., in regime I in Fig. [TJ the equations of state reduce in lowest 
order to 

h-^B = -cf(f c + e)ip A + c£iP A , (36a) 

h-1> A = -cf(/ c + e)V B + cfVl- (36b) 
Solving for ipA and tAb, we obtain 

U=o ~ V>b U=o ~ e, (37a) 

Vu|e=0~V>B|e=0 ~ fr 1/2 . (37b) 

Substitution into Eq. (|35p yields 

J"(e,0) - e 3 , (38a) 

J"(0,/i) ~ /i 3/2 , (38b) 



Percolation on bipartite scale-free networks 



9 



thus, according to Eqs. (f29|) . v c = 3 and Vh = 3/2. All other exponents can now be calculated using 
Eqs. pip . We list the result in the first column of Tabled] 

Next, we focus on the regime in which one of the topological exponents is larger than four, while the 
other exponent has a value between three and four. Without loss of generality, we take in the remainder 
of the text j A as the smallest of the two exponents. The self-consistent equations now reduce to 

h-ipB = -ct(f c + e)ip A + c^Y- 2 , (39a) 
h-^A = -cf(/ c + e)V^ + cfV|, (39b) 



from which we obtain 
1 



1a -3' 



(40a) 



vk = =r^. (40b) 

lA ~ 2 

All other critical exponents are given in the second column of Table [TJ 

Finally, we discuss the exponents in the case 3 < ^ A < 7 B < 4. Although the self-consistent equations 
are slightly different from those in the previous paragraph, we can verify that Eqs. (|40p are still valid. 
The exponents in this regime are thus the same as in the regime discussed in the previous paragraph. 
Therefore, we define region II as the region in which 3 < 7^ < 4 and 7^ < j B , as is illustrated in Fig. [TJ 
The exponents in that regime can be found in the second column of Table [TJ 

4.3.2. Robust networks 

We still have to discuss the exponents for robust bipartite networks. We first focus on regime IV in 
Fig. 1 where both j A and j B are smaller than three. Since f c — 0, the self-consistent equations close to 
the critical point reduce to 

h-ip B = -cf(&p A y A - 2 , (41a) 
h — ipA = -cf(«/>s)~ B ~ 2 , (41b) 

where cf > and cf > 0. These equations can be derived by minimization of 

-=il>Afl>B+ V -hffii - CW i , (42) 

i=A,B ^ ~ 1 

with respect to the order parameters. At the saddle point, we may write: 

f=E-^-|fef^r' (43) 

i=a,b ^ '* ' 

By solving Eqs. (|4"T) for Vi, one obtains 

i> A \h=o ~ e i-(-A-2)(- B -2) t (44a) 

(~ B -l)(- A -2) 

ip B \h=o ~ e i-ciA-2m B -2) f (44b) 

tp A \e^o ~ \h\^/e, (44c) 

V>sUo ~ \h\^/e. (44d) 



Substitution into Eq. (|4"3")l yields 

(7a~1)(7b-1) 
l-(7 A -2)(7 B -2)' 



(45a) 



v h = (45b) 
7s - 2 

The other exponents are listed in the last column of Table 1. 
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Table 1 

Critical exponents in the different regimes for which 7^ < j B . Note that the exponents for the primed 
regimes in Fig. [T] can be obtained by interchanging j A and j B in the expressions for the regimes without 
primes. 





Regime I 


Regime II 


Regime III 


Regime IV 




Fragile 


Fragile 


3 -7 a 
1-7a 


1 ( 3-7a I 3-7 B \ 

2 \ i-iA 1 i-7b ; 




1 


1 


7a-1 


7a-1 




7.4-3 


3-7a 


1-(7a-2)(7b-2) 


T 


5/2 


27,1-3 
7.4-2 


3 


27 B -3 

7b -2 


a 


1/2 


7.4-3 


3 -7 a 


1-(7a"2)(7b-2) 


7A-2 


7.4-1 


(7a-1)(7b"2) 


a 


-1 


5-7 A 


4(7a"2) 


7+37a7b-5(7a+7b) 


7a-3 


3-7a 


1-(7a-2)(7b"2) 


~,> 


1 


1 





(3-7b)(7a-1) 
1-(7a-2)(7b"2) 


8 


2 


7a~2 


1 


7b-2 



Finally, we focus on the regime in which j A < 3 and 7 y B > 3, i.e. regime III on Fig. 1. Although the 
expansion of the self-consistent relations, Eqs. (|T4")) . depends on the actual value of 7 B , the scaling of the 
order parameters does not. After some calculations, we obtain 

,, = (40a) 
u h = 2. (46b) 
We list the other exponents in the third column of Table 1 . 
4.3.3. Discussion 

To summarize, we have calculated the critical exponents of the percolation transition for all physically 
relevant regimes. Note that all exponents in the regimes I, II and IV reduce to the correct expressions 
for monopartite graphs if j A — j B [17] . Note also that, as expected, the usual mean-field results for 
percolation are recovered only in regime I [21|28j . In all other regimes, we find non- universal exponents 
which depend on the constants 7,. 

The only dependence on the bias exponent a arises through the exponents 7 y i as defined in Eq. ((2J. 
Biased percolation with bias exponent a on a bipartite scale-free network with topological exponents 7a 
and 7b thus has the same critical behavior as random percolation (a = 0) on a network with topological 
exponents 7^ and j B . We thus conclude that the results for bipartite graphs are a generalization of those 
for monopartite graphs. 



5. CONCLUSION 

As an extension of previous work on scale- free graphs with a single type of nodes [15117) . this article 
studies the biased removal of links in a bipartite network. A bipartite network consists of two types of 
nodes, with links only running between nodes of unlike type. We assume that the degrees of both nodes 
are distributed according to a decreasing power law, without any correlations or mixing patterns between 
the two types of nodes. This model can be used to describe real-life social networks, such as for instance 
the network of hetero-sexual contacts [13] . In our percolation process, we reinstall a fraction / of the 
links in a single sweep. Moreover, the process is biased in the sense that the link removal probability 
depends on the degrees of the nodes they connect. We attach a weight (kq)~~ a to a link between nodes 
with degrees k and q and retain links with a probability proportional to their weights. In the study at 
hand, the bias exponent a is a positive number smaller than one. Therefore, links between hubs are more 
prone to failure, thus the process is centrally biased. 
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The most important conclusion of the article is that almost all results for the bipartite graphs are 
extensions of the results for monopartite graphs. The robustness or fragility of a bipartite network 
simply depends on the fragility or robustness of the underlying monopartite graphs. Moreover, as for 
monopartite graphs, biased percolation can be mapped to the universality class of random removal on a 
different scale-free bipartite graph. 

To obtain our results, we extended the generating functions theory for monopartite graphs to incorpo- 
rate two types of nodes. Self-consistent equations for the generating functions allowed the determination 
of the critical fraction of the percolation process. Our results show that the critical fraction of a bipartite 
network can be written as a function of the critical fractions of the associated monopartite graphs. If 
at least one of those monopartite graphs is robust, also the bipartite network will be robust. Based on 
these conclusions, we constructed a phase diagram in Fig. 1 picturing the different percolation regimes. 
Also the scaling of the critical threshold as a function of the network size was found to be an extension 
of the scaling of the critical fractions of the monopartite graphs. To calculate the critical exponents close 
to the percolation transition, a finite-size scaling theory was developed using the analogy between edge 
percolation and the q — > 1 limit of the g-state Potts model. The scaling theory was solved by expanding 
the self-consistent equations close to the percolation threshold. Results for the exponents are given in 
Table 1. All exponents only depend on the constant a through the constants j A and j B , defined as 
7 ; = (7^ — a)/(l — a). The critical behavior of biased percolation with bias exponent a on a bipartite 
net with topological exponents 7$ with i = A, B thus is the same as the critical behavior of random 
removal on a network with exponents 7^. Therefore, we extended the main result of the study concerning 
monopartite graphs to bipartite graphs. 

The generating functions theory we introduced can in principle be extended to include general multi- 
partite networks with an arbitrary number of types of nodes. In such networks, nodes can share edges 
with different types of nodes. For each type of nodes, a new generating function must be introduced, 
thereby increasing the complexity of the mathematics greatly. Since every new generating function re- 
quires an additional self-consistent equation, it becomes impossible to extract specific results, except in 
certain limiting cases. Much progress on this scheme has already been worked out in Ref. [12], but a 
calculation of critical exponents still remains an open issue. Further, in Ref. |36j the authors performed 
a study of the extremum events of scale- free networks, thereby focussing on the statistics of the extreme 
connectivities. In this context, it would also be interesting to investigate the evolution of such distribu- 
tion functions during the process of our network reconstruction and more especially near the point of 
percolation. 
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